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ABSTRACT 

m 

y—i Motivated by the need for rapid localisation of newly discovered faint millisec- 

^ ond pulsars (MSPs) we have developed a coherently dedispersed gating correlator. 

o This gating correlator accounts for the orbital motions of MSPs in binaries while 

^ folding the visibilities with best-fit topocentric rotational model derived from 

^ periodicity search in simultaneously generated beamformer output. Unique ap- 
pheations of the gating correlator for sensitive interferometric studies of MSPs are 

^ illustrated using the Giant Metrewave Radio Telescope (GMRT) interferometric 

HH array. We could unambiguously locahse five newly discovered Fermi MSPs in the 

4^ on-off gated image plane with an accuracy of ±1". Immediate knowledge of such 

I precise position allows the use of sensitive coherent beams of array telescopes for 

^ follow-up timing observations, which substantially reduces the use of telescope 

^ time (~ 20 X for the GMRT). In addition, precise a-priori astrometric position 

'— ' reduces the effect of large covariances in timing fit (with discovery position, pul- 

^ sar period derivative and unknown binary model), which in-turn accelerates the 

^ convergence to initial timing model. For example, while fitting with precise a- 
priori position (±1")) timing model converges in about 100 days, accounting the 

^ effect of covariance between position and pulsar period derivative. Moreover, 

^ such accurate positions allows for rapid identification of pulsar counterpart at 

^ other wave-bands. We also report a new methodology of in-beam phase cali- 

^ bration using the on- off' gated image of the target pulsar, which provides the 

>■ optimal sensitivity of the coherent array removing the possible temporal and 

^ spacial decoherences. 

Subject headings: pulsars: general — pulsars: individual (Jl 120—3618, J1207— 5050, 
J1551-0658, J1646-2142, J1828+0625) - techniques: interferometric 
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Introduction 



The ongoing sensitive surveys (e.g. Fermi directed radio searches, Ray et aL (2012); 



Green Bank drift scan Surve}|^ Green Bank North Celestial Cap (GBNCC) survejj^ High 



Time Resolution Universe (HTRU) survey at Parkes, Keith et al. (2010); Pulsar Arecibo 



L-Band Feed Array (PALFA) survey, Cordes et al. (2006)) have discovered a number of 



intriguing fainter millisecond pulsars (MSPs). As a result in last three years the population 
of Galactic disk MSPs is increased by about 45/^ Sensitive follow-up studies of these newly 
discovered MSPs using coherent beams of array telescopes, or at higher frequencies using 
single dishes, are hindered by the large uncertainties associated with the discovery positions. 
For example, discovery uncertainties range from 40' for the GMRT at 322 MHz (GMRT— 322) 
to 4' for the Arecibo at L-band. Sensitive coherent array follow-up observations significantly 
reduces the use of array telescope time (~ 20 x for the GMRT), which is important as 
arrays are the future for large radio telescopes. Such coherent array observations improves 
the uncertainties in time-of-arrivals (TOAs) and allows to generate much closely spaced 
TOAs in order to avoid ambiguous phase connection while timing the binaries with shorter 
orbits. Traditionally long-term pulsar timing programs are used to reduce such positional 
uncertainties, requiring significant amount of observing time. Simultaneous timing fit with 
discovery position and unknown binary parameters can be affected from large covariances, 
specially for long period binaries. In addition, the covariance between position and pulsar 
period derivative (P) limits the convergence of timing fit even with known binary model. 
The effects of such large covariances in timing fit can be minimised with precise a-priori 
astrometric position. 

Being compact objects pulsars are effectively seen as point sources in interferometric 
imaging. Pulsars specially MSPs are weak radio sources in continuum image plane, having 
fluxes in the range of few mJy even at lower frequencies. Identification of pulsar counterpart 
in the continuum image can also be confused with the other sources in the field of view. 
Considering the narrow duty cycle, 3%-10% ( |Henry fc Paik |1969D , the detect ion significance 
of a pulsating point source can be largely improved by removing the off-pulse noise. This is 
achieved through pulsar gating, where the continuum image is sampled synchronously with 
the pulsed signal to generate a number of gated images for different pulse phases (e.g. pulsar 



gating at the ATCA, Lazendic (1999)). Background sky subtracted on-Oj^^gated image allows 



to unambiguously identify the location of pulsed emission (Camilo et al. 2000). However 



■^httpi/Zw-ww. as. wvu.edu/ pulsar/GBTdrift350/ 
"^http: / / arcc.phys.utb.edu/gbncc/ 

^http://astro. phys.wvu.edu/GalacticMSPs/GalacticMSPs.txt 
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such precise position determination using gating is hitherto not been reported for MSPs. 
High time resolution gating requirements for MSPs, with the gating window ~ 100 fis, could 
be computationally challenging. More importantly for MSPs, since period and intrinsic pulse 
width are quite small, dispersion correction is required to be done (unlike normal pulsars) 
before such high time resolution gating, in order to account for large dispersion delay. For 
example, considering a MSP with a dispersion measure (DM) of 30 pc/cc, the dispersion 
delay at 322 MHz GMRT frequency across 32 MHz band is ~ 247 ms, which is many times 
of MSP period. In earlier studies gating has been done after correlation using incoherent 
dedispersion. But incoherent dedispersion smearing across frequency channels can be 5 times 
larger than intrinsic pulse width (considering a 2 ms MSP at 30 pc/cc DM with 5% duty cycle 
for GMRT— 322), significantly reducing the number of effective gates due to pulse smearing. 



Thus reconstructing intrinsic pulse width with coherent dedispersion (Hankins & Rickett 



1975 ) will be beneficial while performing gating with larger bin numbers requiring sufficient 
time resolution. In addition, since majority of MSPs are in binary system, the effect of 
orbital motion on pulsar period has to be accounted for in MSP gating correlator. 

Alternatively localisation of newly discovered pulsars can be achieved by forming mul- 
tiple coherent beams covering the primary beam of the telescope. For example the grating 
response of a linear array like the WSRT can be used for localisation with an accuracy of 



few arcminutes ( Rubio-Herrera et al. 2012). Whereas with continuum imaging followed by 



the multi-pixel beamformer using nonlinear array like the GMRT ( |Roy et al. |[20l2| we have 



achieved the positional accuracy of few arcseconds (determined by the synthesized beam of 
the array). But this method is not efficient to localise relatively fainter MSPs having very 
low detection significance in continuum image plane. 

In order to improve the detection significance in image plane and to achieve positional 
accuracy of the order of 1", we have developed a coherently dedispersed MSP gating correla- 
tor at the GMRT. In addition to aperture synthesis at moderate time constants (~ seconds) 
and high time resolution incoherent and coherent beam formation (~ 30 /is), the GMRT 



Software Back-end (GSB; Roy et al. (2010)) is equipped to stream the raw base-band sam- 
ples from all the antennas to array of storage disks. In this paper we describe design and 
implementation of the MSP gating correlator using raw Nyquist sampled base-band data 
and its application to obtain the precise locations of five newly discovered MSPs from the 
Fermi directed searches. In addition we demonstrate an unique application of gated imaging, 
in using pulsar as a phase calibrator to achieve the optimal sensitivity for the coherent array. 
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2. MSP gating correlator 

Design of tlie MSP gating correlator that can increase the signal-to-noise (S/N) of a 
pulsar in image plane by a factor of 3 to 6 (approximately proportional to the inverse square 
root of duty cycle) is described below. 

(a) Coherent dedispersion : Since baseline based incoherent dedispersion performed 
before folding the visibility time-series at the best-fit topocentric model is not adequate for 
the study of MSPs at lower frequencies, we implement antenna based coherent dedispersion. 
Antenna based coherent dedispersion (~ N operations) can also be computationally favorable 
compared to baseline based incoherent dedispersion (~ operations) for the future arrays 
with large number of elements (N). Coherent dedispersion module is running on the recorded 
raw base-band data prior to correlation, and is parallelized over telescopes. The dedispersed 
voltage samples are written to a shared memory ring buffer. 

(b) Correlation and folding : Dispersion corrected raw voltage samples read concurrently 
from the shared memory are correlated to generate the high time resolution visibility time- 
series. These visibilities are fed into a gating module which bins the data in multiple gates. 
For a given MSP, while folding using best-fit topocentric rotational model, number of gates 
and intermediate time resolution of dedispersed visibilities, are empirically adjusted to obtain 
optimal S/N. In order to retain this optimal S/N in a 1 hr gated image for a 2 ms MSP 
(considering gate width approximately equal to pulse width), fractional period error per 
rotation (P) is required to be below ~ 10~^^, i.e. MSP period (in ms) is required to be 
correct up to 7 significant digits. Folding with a fixed topocentric period is sufficient for 
imaging isolated or loosely coupled binary MSPs. However, we have parameterised the 
rotational model to include period (P), acceleration (P) and jerk (P) for folding MSPs in 
tighter orbits. For newly discovered MSPs this topocentric rotational model is derived from 



PRESTO (Ransom et al. 2002) based periodicity search for the same observation using 
simultaneously generated incoherent beamformer output. Whereas for MSPs with known 
ephemeris we have used TEMPO|^ based predictions. The folded visibilities on each gate are 
then integrated up to a time resolution of 16 seconds and the complex correlation output 
from each gate are finally written to disks for further processing. 

(c) on-Oj(f gating : The on-pulse gate is also identified using the coherent, or incoherent 
where the position of the MSP is poorly known, beamformer output, as the arrival time of 
the pulse is unknown a-priori for new pulsars. The gates containing the pulsed emission are 
grouped to form the on-pulse data, and the off-pulse data is formed with same number of 



^littp://www.atnf .csiro.au/research/pulsar/tcmpo 
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gates covering the off-pulse phases. In order to unambiguously identify location of pulsed 
emission we have generated on-off visibility bin. In addition, by subtracting the nearby 
off-gate from the on-gate underlying systematics generated from instrumental effects as well 
as from radio frequency interferences are canceled out, resulting in further improvement in 
noise statistics. 

(d) Calibration and imaging : Coherently dedispersed folded on- Oj^ visibility data are 



then flagged for removing the outliers and calibrated for solving the complex gains ( Prasad & 
Chengalur||2011 ). Calibrated visibilities are imaged and deconvolved using standard imaging 
package AIPS[^ MSP is the only source in this on-ojf gated image, considering that on and 
off gates contain equal flux for other sources. 



3. Localisation of MSPs using gating correlator 

MSP gating correlator is efficient in localising fainter MSPs, where fluxes are around 1 
mJy resulting in very low detection significance on continuum image plane (rms noise in 1 hr 
at GMRT— 322 ~ 500 /iJy). Positional accuracy achieved from on- ojf gaied image depends 
upon hour angle of the observations and S/N of pulsar detection. Considering GMRT— 322 
synthesized beam (FWHM) ~ 10", the positional accuracy of a MSP scales according to 
FWHM/(2xS/N). For a typical S/N of 5 (Table. |2|, an accuracy ±1", is determined by 
the AIPS task JMFIT. Such a-priori astrometric accuracy accelerates the convergence in 
pulsar timing. The newly discovered MSPs have large uncertainties in the a-priori position, 
requiring timing span of the order of a year to overcome the effect of covariance between 
position and P. However, while fitting with more precise a-priori position, having ±1" 
uncertainty, the positional accuracy of ~ 1 mas and a convergence in detection of P (with 
AP/P ~ 0.04) can be achieved in only about 100 days. 

We performed gated imaging for five MSPs discovered in Fermi directed radio searches. 
Among those, PSR J1120-3618, J1646-2142 and J1828+0625 were discovered by us in 



GMRT-322, whereas PSR J1207-5050 in GMRT-607 (Ray et al. 2012) and PSR J1551-065 



in GET— 350 (Bangale et al. 2012[ ). The parameters of these MSPs such as period, DM 



and mean flux are listed in Table [TJ Quoted mean flux is obtained using the simultaneous 
incoherent beamformer output. For PSR J1207— 5050 the gating observations were done at 
607 MHz while rest of the MSPs were observed at 322 MHz. 

PSR J1120-3618 is a serendipitously discovered MSP in the GMRT-322 field-of-view 



^http: / / www.aips.nrao.edu 
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of known Fermi MSP Jl 124— 3653. This is a very faint MSP (mean flux ~ 300 yuJy) in a 
relatively tighter orbit. PRESTO based search pipeline reports an acceleration equal to 0.1 
m/s^, corresponding to a significant measurement of P equal to 2xl0~^^ s/s. This introduces 
a period (in ms) error at 1 part in 10^ over 3 hrs of observation, whereas for optimal folding 
an accuracy up to 1 part in 10^ over 3 hrs is needed. Thus accounting for acceleration 
is required to retain the S/N of this MSP in on-Oj^ gated image plane. The pulse phase is 
binned in 11 gates (gating window ~ pulse width) at the intermediate time resolution (491.52 
/xs) of the dedispersed visibilities. The on-Oj(f gated image for this MSP is shown in Fig. [l] 
Interestingly this MSP is located at 57' offset from the pointing center, which is out-side the 
GMRT— 322 beam-width (~ ± 40'). To localise the MSP we have performed multi-faceted 
dPerley |p99| gated imaging which is also corrected for primary beam effect as the pulsar is 



located at the edge of the beam. In order to obtain the gated image of the full field-of-view, 
separate facet images are interpolated and averaged onto a larger grid using the AIPS task 
FLATN. A lO'xlO' facet of the on-Oj^ gated image shows the MSP (Fig. [T]) with 5a detection 
significance, resulting in gated flux of 1.2 mJy which is 4x the mean flux as expected from 
gating. The parameters related to gated imaging are listed in Table. [2j Sensitive coherent 
beam is formed towards the pulsar location by steering the phase center and an expected 



sensitivity improvement ~ 4x of incoherent array detection is achieved (Roy et al. 2012). 



PSR J1207-5050, J1551-0658, J1646-2142 and J1828+0625 are also successfully lo- 
calized in their respective on-Oj^ gated images (shown in Fig. [2]). Table. |2] lists the related 
parameters. PSR J1207— 5050 and J1646— 2142 are found within the error radius of the as- 
sociated Fermi sources. Whereas PSR J1551— 0658 and J1828+0625 are located at an offset 
of 20' and 26' respectively from their pointing centers (outside the error radius), indicating 
that these MSPs are not likely to be associated with the corresponding gamma-ray sources. 
Large primary beam of the GBT— 350 and large incoherent beam of the GMRT— 322 have 
allowed such serendipitous discoveries. The gated fluxes for all these four MSPs are 3 to 6 
times of the respective mean fluxes as expected from the S/N improvement by gating. 



4. In-beam phase calibration using pulsar 

The sensitive coherent beamformer allows study of high time resolution temporal vari- 
ations of pulsars. In order to form coherent beam using an array telescope like the GMRT, 
antenna based complex gains (amplitudes and phases) need to be solved using recorded 
visibilities on a calibrator source. The optimal baseline length over which the array can 
be coherently added is limited by the perturbations in ionospheric phases, which are more 
severe at lower frequencies. Coherent array sensitivity degrades with time due to the tem- 
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Table 1: Parameters of the concerned MSPs 



PSR 


Period 


Dispersion 


Mean 




(ms) 


measure (pc/cc) 


flux" (mJy) 


J1120-3618 


5.55 


45.1 


0.3 


J1207-5050 


4.84 


50.7 


0.2 


J1551-0658 


7.09 


21.6 


1.0 


J1646-2142 


5.85 


29.7 


2.1 


J1828+0625 


3.63 


22.4 


1.3 



a : Mean flux of PSR J1207-5050 is measured at 607 MHz, rest are in 322 MHz 



Table 2: Parameters related to gated imaging 



PSR 


Gated J2000 


Offset from 


Number 


Observing 


Gated 


Gated 




position 


pointing 


of gates 


duration 


flux" 


SNR 




(Errors in ") 


centre 




(min) 


(mJy) 




J1120-3618 


11^20°^22!405 (I'.'l); 


57' 


11 


180 


1.2 


5 




-36°18'32'.'17 (2'.'2) 












J1207-5050 


12'^07"'21!811 (0'.'4); 


6.2' 


10 


120 


1.1 


6 




-50°50'30'.'27 (1'.'4) 












J1551-0658 


15^51'^07!215 (C.'a); 


20' 


14 


60 


5.8 


11 




-06°58'06'.'51 (0'.'6) 












J1646-2142 


16M6'^18n27 (0'.'9); 


10' 


12 


60 


11.3 


11 




-21°42'08'.'96 (1'.'4) 












J1828+0625 


18^^28^28^030 (I'.'O); 


26' 


15 


45 


3.5 


6 




06°25'00'.'52 (1'.'3) 













a : Gated flux of PSR J1207-5050 is measured at 607 MHz, rest are in 322 MHz 
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Fig. 1. — Localization of the PSR J1120— 3618 using MSP gating correlator. The MSP 
(marked in the image) is detected in the on-off gated image at 57' offset from pointing 
centre. 



poral decoherences caused by instrumental as well as ionospheric phase fluctuations. Such 
degradation can be reduced with interleaved calibrator observations. However, applying the 
phases derived from a distant calibrator into a pulsar field can cause some more decoherences 



due to underlying different ionospheric inhomogeneities (Thompson et al. 1986). We have 



used the on-Oj^ gated image of the target pulsar as a sky model to solve for antenna based 
residual stochastic phase errors (affecting the data on short time scale) as well as the broad 
band phase offsets applying self-calibration in AIPS. This process produces a set of phase 
solutions with time, written in SN table generated by AIPS. While forming the coherent 
array, these residual phase solutions are recursively applied, in addition to initial narrow- 
band phase corrections derived from the calibrator visibilities. Background sky subtracted 
on-Oj(f gated image of the target pulsar provides a better model to solve for phases since the 
effects of instrumental as well as RFI artifacts are canceled out. A pulsar with about 10 
mJy flux density allows this in-beam calibration process to converge (ensuring more than 
3cr detection in each frequency channel during calibration with 10 minutes cadence). This 
in-beam calibration is the optimal way of coherently adding all the working antennas of an 
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Fig. 2.— The on-ojf gated images for PSR J1207-5050, J1551-0658, J1646-2142 and 
J1828-I-0625. All the MSPs are marked in the respective lO'xlO' facet images. 



array telescope. 

We have generated on-Oj^ gated image for PSR B1804— 08 using 1 minute of base-band 
data, which is shown in upper left panel of Fig. [3j With the in-beam calibration a sensitivity 
improvement of 3.5 x compared to the conventional coherent array is seen in the folded profile 
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(upper right panel of Fig. [s]). This improvement includes contributions from the increased 
coherence length of the array (achieved from 50% increase in number of antennas) and better 
modeling of phase errors. Since the spectral voltages from the antennas are optimally added 
in phase, the sensitivity improvement can also be visualized as a reduction in the spectral 
noise in the dispersed pulse phase versus frequency plot (bottom panel of Fig. |3]) generated 
with DSPSR (Straten & Bailes 2010). In addition, the temporal decoherence affecting the 
long observing scans of target pulsars can be avoided by recursive in-beam calibrations at 
short time scale, without slewing to a distant calibrator location. Use of pulsar as an in-beam 
phase calibrator can introduce a lateral shift in the image domain while pulsar is not at the 
phase centre. This can be calibrated using known position of any in-field catalogue source 
allowing to perform an astrometric measurement of pulsar. 



5. Summary and future scope 

We report design and implementation of a coherently dedispersed MSP gating correlator 
which accounts for a full timing model including orbital motion while folding the visibility 
time-series. In this paper we have unambiguously determined the precise positions, with ±1" 
accuracy, for five newly discovered MSPs using the on- off gated imaging. Locatisations of 
such relatively fainter MSPs are greatly benefited from the significant enhancement in S/N 
on the gated image plane compared to the normal synthesis observations. Knowledge of such 
accurate positions in immediate observations after discovery allows follow-up observations 
with sensitive coherent array, substantially reducing the use of telescope time by an order 
of magnitude. Inaccuracy of the astrometric model, associated with large uncertainties 
in discovery positions, increases the length of data span required to reduce the effect of 
covariance between position and P in timing fit. While determining an unknown binary 
model, large covariance between position and the binary parameters, specially for long period 
binaries, can influence the timing flt, which can be minimised with accurate astrometric 
position (~ ±1"). In addition, for pulsars located near the ecliptic plane, the astrometric 



inaccuracy in ecliptic latitude from the timing flt is much larger (Lorimer & Kramer 2004) 
indicating the need for interferometric gated observations. This accurate localisation will 
also facilitate the search for pulsar counterpart and possible binary companions at optical 
and X-ray. Current astrometric positional accuracy (±1") is decided by array size and S/N 



of detection. Application of MSP gating correlator in the VLBI (Deller et al. 2009) and 
the SKA scale improves this a-priori astrometric accuracy (e.g. at 1.4 GHz, 3000 km SKA 
baseline can give astrometric precision as 15 /xas ( Smits et al. | 2011)). 



In addition we have used the on- off gated image of the target pulsar to correct the 
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Fig. 3. — In-beam phase calibration using PSR B1804— 08. Upper left panel shows the on-off 
gated image with a zoom around the pulsar. Upper right panel shows the coherent array 
sensitivity improvement in the folded pulse profile using in-beam phase calibration. The 
dispersed pulse phase with frequency is shown in the bottom panel. The similar sensitivity 
improvement with reduction of spectral noise (S/N are mentioned on the right corners of the 
respective plots) is seen on the right panel with respect to the left panel. 



residual phase errors in order to avoid degradation of coherent array sensitivity caused by 
decoherences from instrumental and ionospheric phase fluctuations. Such in-beam phase 
calibration using a target pulsar ensures optimal sensitivity of the coherent array. We believe 
that this methodology will be very fruitful for recently commissioned and up-coming array 
telescopes (e.g. LOFAR, |de Vos et al~~| ([2009|; MWA, [Lonsdale et "aT"! (12009]); ASKAP, 



Johnston et al. (2008); MeerKAT, Jonas (2009)). 



-VI- 



Even though the design is primarily motivated by the requirement of locahsation of 
newly discovered Fermi MSPs, this gated imaging has the potential to unfold some other 
interesting properties of MSPs. Firstly, the MSP gating correlator will allow to perform 
independent measurement of parallaxes and proper motions even for relatively fainter binary 
MSPs, which will largely benefit pulsar timing and will probe the interstellar medium (ISM) 



at various line-of-sight (McGary et al. (2001), Lorimer & Kramer (2004), Smits et al. 



I ( 2011[ )). Secondly, the study of un-pulsed emission associated with pulsar wind nebulae 



(PWN), can be performed by imaging the pulsar field when pulsed emission is off (Gaensler 



et al. I 2000). Angular extent of PWN for high energy pulsars in low density ISM can 
be as large as few arcminutes (e.g. Frail et al. (1994)), which can be detected with low 



frequency gated imaging using this MSP gating correlator. Finally, pulsar itself may have a 
weak off-pulse emission, having magnetospheric origin, coming from close to light cylinder 



(Perry a& Lyne (1985), Basu et al. (2012)). Origin of possible off-pulse emission from 
MSPs can be probed with the gated images as a function of pulse longitude and any possible 
co-location with the gamma-ray emission region will be interesting to investigate. Thus the 
design of coherently dedispersed MSP gated imaging assures broader scientific returns, while 
also having importance for SKA applications. 



6. Acknowledgments 

We thank computer group at GMRT and NCRA. We thank Dr. David Thompson and 
Dr. Paul Ray for helping us to join PSC and Prof. Yashwant Gupta and Prof. Dipankar 
Bhattacharya for being members of the GMRT team. We also thank Prof. Ue-Li Pen for 
insightful discussions and Dr. Maura Mclaughlin for inspiring to work on J1551— 0658. We 
thank the anonymous referee for comments that improved the quality of the paper. We 
acknowledge support of telescope operators. The GMRT is run by the National Centre for 
Radio Astrophysics of the Tata Institute of Fundamental Research. 

REFERENCES 

Bangale, P. et al, 2012 (in preparation). 

Basu, R., Mitra, D., & Athreya, R., 2012, ApJ, 758, 91. 

Camilo, F., Kaspi, V. M., Lyne, A. G., et al, 2000, ApJ, 541, 367. 

Cordes, J. M., Freire, P. C. C., Lorimer, D. R., et al. 2006, ApJ, 637, 446. 



Deller, A. T., Tingay, S. J. & Briskan, W., 2009, ApJ, 690,198. 

de Vos, M., Gunst, A. W., & Nijboer, R., 2009, IEEE Proc, 97, 1431. 

Frail, D. A., Goss, W. M.,& Whiteoak, J. B. Z. 1994, ApJ, 437, 781. 

Gaensler, B. M., Stappers, B. W., Frail, D. A., Moffett, D. A., & Johnston, S., 2000, MNRAS, 
318, 58. 

Hankins, T. H. & Rickett, B. J, 1975, Pulsar signal processing. Methods in computational 
physics, 75, 55. 

Henry, G. R., & Paik, H. J., 1969, Nature, 224, 1188. 

Johnston, S.; Taylor, R.; Bailes, M., et al., 2008, Experimental Astronomy, 22, 151. 
Jonas, J. L., 2009, Proceedings of the IEEE, 97, 1522. 

Keith, M. J., Jameson, A., van Straten, W., et al., 2010, MNRAS, 409, 619. 

Lazendic, J. S., 1999, Australia Telescope Compact Array Users Guide (Sydney : ATNF). 

Lonsdale, C. J., Cappallo, R. J., Morales, M. F., et al. 2009, IEEE Proc, 97, 1497. 

Lorimer, D. R., & Kramer, M., 2004, Handbook of Pulsar Astronomy, Vol. 4. Cambridge, 
UK, 211. 

McGary, R. S., Brisken, W. F., Fruchter, A. S., Gross, W. M., & Thorsett, S. E., 2001, ApJ, 
121, 1192. 

Perley, R. A., 1999, ASP Conference Series, 180, 383. 

Perry, T. E., & Lyne, A. G. 1985, MNRAS, 212, 489. 

Prasad, J., & Chengalur, J. 2011, Experimental Astronomy, 141. 

Ransom, S. M., Eikenberry, S. S., & Middleditch, J. 2002, AJ, 124, 1788. 

Ray, P. S. et al., 2012, 2011 Fermi Symposium proceedings, |arXiv: 1 205. 3089} 

Roy J., Gupta Y., Ue-Li Pen, et al., 2010, Experimental Astronomy, 28, 55. 

Roy J., Bhattacharyya B., & Gupta Y., 2012, MNRAS Letters, 27, 90. 

Rubio-Herrera, E., Stappers, B. W., Hessels, J. W. T., & Braun, R., 2012, MNRAS,tmp272. 



- 14 - 



Ryle, M., & Vonberg, D., 1946, Nature, 158, 339. 

Smits, R., Tingay, S. J., Wex, N., Kramer, M., & Stappers, B., 2011, A&A, 528, 108. 
Straten, W. van, & Bailes, M., 2010, [arXiv: 10 08739731 

Swarup, G., Ananthakrishnan, S., Subrahmanya, C, R., et al., 1997, in High Sensitivity 
Radio Astronomy, ed. Ja ckson, N., and Davis, R., J., (Cambridge: Cambridge Uni- 
versity Press), 217. 

Thompson, A.R., Moran, J.M., & Swenson, G.W., 1986, Interferometry and Synthesis in 
Radio Astronomy, John wiley and Sons, New York. 



This preprint was prepared with the AAS I^TjTpC macros v5.0. 



